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Motivated by the tentative observation of superluminal neutrinos by the OPERA experiment, we 
present a model of active-sterile neutrino oscillations in which sterile neutrinos are superluminal and 
active neutrinos appear superluminal by virtue of neutrino mixing. The model demonstrates some 
interesting possibilities and challenges that apply to a large class of models aiming to explain the 
OPERA result. 

I. INTRODUCTION 

<N 

O ■ Measurements of the arrival times of muon neutrinos from the CERN CNGS beam at the OPERA detector 730 km 
away suggest that they travel superluminally with (v — cY/c = (2.37 ± 0.32(stat.) ± 0.29(sys.)) x 10 -5 Many 
U J • interpretations of this stunning result have been proposed [2| . 

' Inspired by the OPERA anomaly, we present a specific realization of a class of models that may be viewed either as 
superluminal travel of a gauge-singlet sterile neutrino via extra-dimensional shortcuts [3| or alternatively as Lorentz 
violation for sterile neutrinos as viewed from our four-dimensional spacetime Q. We emphasize at the outset that it 
is not our intention to explain the OPERA data, but to simply provide a concrete model of superluminal neutrinos. 
As described at length in Ref. Q , a superluminal sterile neutrino is well-motivated within the context of brane- 
Oh| world phenomenology. The active neutrinos, carrying electroweak gauge charge like all other Standard Model (SM) 
fields, are described as open string excitations with their string endpoints confined to our 3 + 1 dimensional brane. 
On the other hand, the sterile neutrino, carrying no gauge charge, is characterised by a closed string, free to roam 
the extra-dimensional bulk as well as the brane (in the fashion of the "gauge-singlet" graviton). Thus, its geodesic 
between two points on the brane will include travel in the bulk. The net result in general will be a shorter transit 
distance; such shortcuts were proposed a decade ago for gravitons in Ref. From the point of view of our brane, 
the sterile neutrino will appear to travel superluminally. An analogy would be a comparison of the light transit time 
ly-j , and distance when confined within a curved optical fiber, and the light transit time and distance when traveling 
the straight path between the fiber endpoints. As proposed in Ref. the shorter distance through the bulk could 
be a result of brane fluctuations within the bulk. These fluctuations could be thermal, gravitational, or quantum 
mechanical in origin. Also, the difference in the limiting velocities between active and sterile neutrinos 5v is related to 
the geometry of the brane fluctuation. The relation Sv — (4^) 2 was found, where A is the (classical) amplitude of the 
brane fluctuation in the bulk direction, and k is the wave number of the brane fluctuation along the brane direction. 
Thus, Sv is basically the dimensionless aspect ratio of the brane fluctuation. 
■ i-h , This article is organized as follows. We derive the oscillation probabilities when only one active-sterile mixing angle 
is nonzero. In doing so, we extend the results of Ref. from two to three active neutrinos, plus one sterile neutrino. 
We then briefly mention aspects of the OPERA data in the context of our model. Finally, we conclude. 
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II. FORMALISM 



The quantum mechanics of the model is simple. The flavor-oscillation amplitude for a propagating neutrino is 

A(y a ^v fi ) = {v fi \e^ m \v a ). (1) 

A component of Ht that is proportional to the identity cannot affect flavor change, and can be subtracted. We write 
the remainder as 8(Ht) = (5H)t + H(8t) under the assumption that it is small. We are left with 

A{v a ^v p ) = {v p \e-W H » +H W\v a ). (2) 
As in standard oscillations, SH is diagonal in the mass-basis, and at lowest order is equal to 



SH = — diag(m?,m§,---)- 



(3) 
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Upon inserting complete sets of mass eigenstates before and after e~' l ^ SH ' t in Eq. ([2]), the first term there becomes 
U*^ Upj e _l "2B"; the usual definition of the bases-mixing matrix, 



or equivalently, 



(4) 



has been employed. 

A nonvanishing value for the second term in Eq. ((2|) is unconventional, and occurs if the propagation times for the 
neutrino states are not universal. Such a theory assigns different "light-cones" to different states, thereby breaking 
Lorentz invariance. Conversely, a large class of models with Lorentz Invariance Violation (LIV) has been shown to be 
phenomenologically equivalent to state-dependent limiting velocities We note that with differing velocities, one 
has 5t = S(L/v) = —L5v/v 2 , which is —L5v to lowest order (in natural units). In the picture where gauge-singlet 
states are closed strings free to roam the bulk, the limiting velocities are assigned to the flavor eigenstates rather 
than to the mass eigenstatesQ Such a choice affects the equivalence between a sterile flavor traveling the shortened 
geodesies available in the bulk, and a Lorentz- violating, superluminal limiting velocity for the sterile state as viewed 
from the brane. The second term in ([2]) as written is already in a diagonal basis, and 



St = diag(<5i Q , 8t@, ■ ■ ■ ) = — L diag(5u Q , dvp, ■ ■ ■ ) . 



(5) 



It is conventional to put the physics into a Hamiltonian framework. The effective neutrino Hamiltonian in the flavor 
basis is 



Hi 



(F) 



( 


m\ 
















m\ 








\ 








m\ 














rn 



U ] -E 



Svi 














Sv 2 














Sv 3 














5«4 



(6) 



In general, the 4x4 mixing matrix U consists of six angles (the number of planes in four dimensions) and four phases. 
To simplify the analysis, we neglect the three new phases, and for now, set to zero the rotation angles in the 4 — 2 
and 4—1 planes. By keeping the #34 angle in R34 nonzero, we retain the basic features of the model. We have 
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V 
1 



1 
x 1 1 

R 34 



(7) 



in the absence of the new term proportional to Sv's. Here, V is the usual PMNS mixing-matrix among the three 
active-flavor neutrinos, and 



R 



cos 6*34 sin #34 
— sin #34 cos #34 



(8) 



We next write A 



il — m§, and neglect the light masses m 2 pj — 1, 2, 3 relative to 771,4. We assume that the active 
neutrino flavors have the usual limiting velocity c, whereas the sterile flavor has a limiting velocity Sv = Sv^ > 0. 
This seems to us to be the most economic and intuitive application of possibly-differing limiting- velocities. The sterile 
state is qualitatively different from active states in that it has no gauge interactions, and therefore is unconstrained 
by gauge symmetries. We provide more discussion of a qualitatively different sterile neutrino below. 
With these assumptions, the effective Hamiltonian in ^ may be written as 
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(9) 



The qualitative features of Ht F \ in Eq. Q provide for an interesting discussion. At sufficiently low energies, the first 
term on the right-hand-side of -ff(F) dominates, and oscillations proceed in the standard way. The second term on the 
right-hand-side, diagonal in the flavor basis, has an analogy with the famous MSW matter-term. At sufficiently high 



An alternative model arises if one assigns the limiting velocities to the mass eigenstates. In such a model, the mass-squared matrix and 
the Sv matrix are diagonal in the same basis, and so there is no brane-bulk resonance arising from diagonalization of the Hamiltonian 
of Eq. (O below. Yet another possibility is to assign limiting velocities to velocity eigenstates. 
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energies, the eigenstates of the Hamiltonian are nearly flavor states, and oscillations are very suppressed. At some 
intermediate value of energy, the two terms are comparable, and resonance enhancement of the mixing angles may 
occur (if the mixing angle can reach the maximal- mixing value of 45°, as discussed below). 
The matrix in brackets in Eq. @ is equal to 
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(10) 



and is diagonalized by the rotation R34 through an angle 6*34 given by 

sin 26*34 



tan 29 = 



cos26> 34 - 2E 2 5v/A ' 



or equivalently, by 



sin 2 29 



sin 2 26»34 



sin 2 26> 34 + (cos 26*34 - 2E 2 5v/A) z 



(11) 
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Because of the second term on the right-hand-side of Eq. @, we are led to a diagonalization matrix of the form given 
in Eqs. ([7]) and ©, but with 6*34 in Eq. ([5]) replaced by 9. Thus, the matrix which diagonalizes the full Hamiltonian 
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i? 34 (S0 



/ V e i Ve2 Ve3 cos 9 V e3 sin 9 \ 

V^i Vju.2 cos 9 VJi3 sin 

V T \ V T 2 V T 3 cos 9 V T 3 sin9 

\ -sin^ cosf? / 



(13) 



Resonant mixing occurs when the two diagonal elements in Eq. (|10p are equal, i.e., when 



A COS 2034 



2Sv 



In terms of E-r, equations (fTl"]) and (fT2j) may be written as 
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tan 26*34 
1-(JL 



(14) 
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and 



r 2 29 = 



sin 2 26* 
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(16) 



The energy- dependent angle 9 is obtained by taking the inverse sine of Eqs. ()T2j) or (|T6|) . or the inverse tangent of 
Eq. (fTTj) or (IT5|) . Care must be taken to ensure that 9 is chosen in the first octant for E < Er, and in the second 
octant for E > E^,. The functions sin 9 and cos 9 are then readily obtained. Since cos 26*34 is positive definite for 
small 6*34, resonance can occur only if A and Sv have the same sign. Cosmological limits on neutrino masses disallow 



3 eV, so A must be positive. Thus, resonance is possible only if 6V4, > 0. One possibility is 



to have limiting velocities W4 = c, Vi < c [i < 4). The other possibility, more natural in the brane-bulk scenario and 
assumed above; is to have Vi = c (i < 4) and v, > c Q. This latter possibility is discussed more below. 

There are two distinct qualitative differences between the LIV resonance inherent in Eq. (j9|), and the MSW matter- 
resonance. Firstly, The LIV term here grows with energy, whereas the matter term in the MSW Hamiltonian does 
not. Consequently, the LIV resonance will be narrower than an MSW resonance. In other words, a measurement of 
the full width at half maximum (FWHM) may be a signature of the LIV resonance. Secondly, the LIV resonance 
does not violate CPT, whereas the MSW resonance necessarily does; the LIV resonance will occur identically in both 
neutrino and antineutrino channels, in contrast to the MSW resonance. 
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The eigenvalues of -HVf) are 



A 



Ai = A 2 = 0, A 4/3 = A± = — 1 - cos2034 -fr- ± 
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and the eigenvalue differences 5H k j = A& — A? are 
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(18) 



With these eigenvalue differences and the mixing matrix U, we have all the ingredients to obtain all possible 
oscillation probabilities Furthermore, in the model aspresented, there are just three parameters beyond the standard 
three-neutrino parameters. These are A, 6*34, and £r 

The general oscillation formulae are 



P( Va -n/ /l ) = ^-4^ m.{Upj U; k U* a3 U ak } sin 2 

j<k 

which on ignoring phases in U becomes 
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P(v a up) = Sa/3 -4)^Upj U{ )k U aj U ak sin 2 
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(21) 
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For the present case we get 



P{y a -> vp) = 5 a p - 4 x < + sin 
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A relevant variable for neutrino oscillations is L/E. Just as £r, sets the energy scale for the resonance, the length 
scale for the n-th maximum at resonance is set by an interplay of the various 

7r(2n - 1) 
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jk\ 



(24) 



Substituting E = £r into Eqs. (jTTJ) and (fT5|) . we find the values of L/E at E = Er for the three contributing 
amplitudes to be 
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(25) 
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2 In the three-neutrino model of Ref. [j| , consisting of two active and one sterile neutrino, the "1" state is absent and so the v\ row and 
column of U is absent, and V is effectively replaced by the 2x2 matrix R2'i{d,). Consequently in Ref. [H, SH34, SH42, and SH32, are 
all of the same magnitude in the resonance region. 

3 More general LIV scenarios lead to dispersion relations of the form 

where Eq denotes some typical energy scale. The case under discussion correponds to n = 1. For arbitrary n, 

~ 1 cos 2#34 \ ™TT 

~26v ) " (20) 

and the corresponding sin 2 28 and SH k j are obtained by replacing (E/E-n) 2 by (E / En) n+1 in Eqs. dl5t -( fT8l) . 
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Note that the 1—2 submatrix of U is the same as that of V. The matrix V, like U, is unitary. Thus, Ylj=i i U<xj Upj — 
da/3 ~ K*3^/33- Making this replacement, and using the explicit matrix entries in the third and fourth columns of 
Eq. (fl3|) . we arrive at simpler expressions for the three relevant cases: active neutrino survival, active-to-active neutrino 
conversion, and active-to-sterile conversion. Denoting the sterile neutrino by v s and active flavors by v a , v^, • • • , the 
active neutrino survival probability is given by 



•4V?« x 




sin 2 cos 2 I 



v a3 



sin 2 9 Q.-V&) 
cos 2 9 (1-V&). 



(27) 



The active-to-(diffcrcnt) active neutrino conversion probability is given by (and mind the minus sign on the first term 
in brackets) 



P{v„ 
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The active-to-sterile conversion probability is given by 

P(v n -> v. 



sin 



sm 




V% sin 2 26» sin 2 



sin 2 cos 2 9 
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(28) 



L(A + -A_; 



(29) 



Note that correct limits are respected here. Far above the resonance, cos 2 9 and A + approach zero (while sin 2 9 
and A- do not). Thus, each term in the above probabilities vanishes far above -Er, and the sterile state effectively 
decouples, as it must. 

The analytic formalism presented here fails if more than one 6^4 is taken to be nonzero, for then the eigenvalues 
must be derived from a matrix larger than the 2x2 subblock given in Eq. (1101) . However, the formalism goes through 
when a 9j4 other than #34 is taken to be nonzero. We have really described three models here, characterized by a 
nonzero #34, #24, or #14. In Eqs. (|2"7| - (f2"9")) . one need only replace the subscript "3" by "2" or "1" to obtain the 024 
and #14 models, respectively. 

For the 6I34 model, we have 4V 2 3 V 2 3 = sm 2 (29 13 ) sin 2 23 , 4V r e 2 3 y T 2 3 = sin 2 (20 13 ) cos 2 23 , and AV 2 3 V 2 3 = 
sin 2 (20 23 ) cos 4 0i3, for the prefactors to P(v e «-» v^), P{y e «-» i/ T ), and P(v^ f-> v T ), respectively. 

The formalism needs to be extended to be relevant for data away from the resonance. Continuing with the simple 
model with just one nonzero Rj4, we see that a factorization occurs between the squared elements of V and R in U, 
viz. (with no sum on j implied) 
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(30) 



The final form results because only the j th active flavor and the sterile state appear in the Rj4 matrix. This result is 
simple matrix multiplication of the matrix with V-squared elements and the matrix with i?-squared elements. It is 
useful to present the squared elements of the V and R matrices. The i?-squared elements are 



/ 1 
'OIO u 

cos 2 sin 2 
V sin 2 cos 2 / 



(31) 



with obvious generalizations for the other choices of nonzero R±j . We next discuss the ^-squared matrix. 

A PMNS matrix consistent with most neutrino data is the tribimaximal matrix Q, which extended to the 4x4 
case is: 



^4x4 — —?= 
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/ 2 V2 \ 
-1 y/2 

-1 V2 -Vz 
V V6 j 



(32) 
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The ^-squared elements of this matrix are 



/ 4 2 

12 3 

12 3 

\ 6 



(33) 



For example, in the ^34-model, f7 2 3 = [V^] 2 [(i?34)33] 2 = \ cos 2 9. This can also be seen from Eq. (fT3"j) . A unitary 
extension of the tribimaximal matrix to nonzero U e3 is given in Ref . [8J . We find that the conversion probability for 
v ii — > to lowest order in |[/ e 3| 2 , is 

£(A V A ' } ) sin 2 9 cos 2 9 
^) sin 2 9 (34) 

^-j cos 2 6. 

Note that the different choices of the nonzero 9j4 result in just a change in the overall magnitude of the flavor-changing 
probability. 

To study the influence of the new mass-squared scale A on long-baseline and atmospheric data, the atmospheric 
mass-squared scale A atm = \m 2 — m\\ must be entered into the probability formulae. This is done simply in the 
following way: In Eq. (J9j> , when #34 is chosen for mixing the diagonal mass matrix is replaced by 
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(35) 



and when #24 or #14 are chosen for mixing, the diagonal mass matrix is replaced by 
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(36) 



These choices ensure that A atm is not mixed into the resonance condition, making the extraction of the four eigenvalues 
simple. We neglect effects of order A atm /A. The two sign choices correspond to normal and inverted hierarchies, 
respectively, for the active mass spectrum. Notice that because one of the three active states is distinguished from the 
other two, due to its mixing with z/4, the two hierarchies yield different physics. This is analogous to the the MSW 
situation in matter. 

For the $34 model, the ordered eigenvalues of the full Hamiltonian are Xj — {=FA atm /2iJ, =FA atm /2£?, A_, A+}; there 
are three distinct nonzero eigenvalue differences SHkj = \k — Xj, as in the A atm = case discussed earlier. To get the 
oscillation probabilities, these eigenvalue differences SHkj = \k — Xj are inserted into Eq. (|22|) . 

For the $24 and 6*14 models, the eigenvalues are Xj — {0 and A_ in the appropriate order, followed by ±A atm /2i5 
and A+}; here there are six distinct eigenvalue differences. After insertion into Eq. (f2"2"|) . one obtains the oscillation 
probabilities after some straightforward but tedious algebra. 



III. OPERA 



OPERA has inferred a mean arrival time for muon neutrinos of ~ 60 ns faster than the theoretical light-travel 
time, over the 730 km pathlength. The short-bunch (3 ns) beam data newly acquired by OPERA provide a stringent 
constraint on models in which neutrino flavor states or mass states travel with different limiting velocities. This is 
because the new data show an arrival time rms of ~ 16 ns about the mean, with no events arriving within 37 ns 
of the theoretical light-travel time. It appears that there is but one neutrino speed in the data. A fit to these new 
data l 9*] requires the fraction of superluminal neutrinos to be at least 80% at a 3a confidence. (Resonant conversion 
offers the possibility of a 50% mixture on average, while adiabatic conversion and back-conversion can theoretically 
attain 100% conversion to the fast species.) This implies that atmospheric neutrinos oscillate primarily into sterile 
neutrinos, a conclusion that is excluded by Super-Kamiokande. 
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Another daunting constraint comes from the L/E distribution of Super-Kamiokande's atmospheric muon neutrino 
events [l(J. A large scale (compared to 2.5 x 1CU 3 cV 2 ) in the difference of eigenvalues of the Hamiltonian, as 
occurs in active-plus-sterile neutrino models (cf. Eqs. (fTT)l . (fT8]) and (|2T)l ) would provide an energy- independent, 
averaged contribution for the L/E distribution, in violation of these data, unless the active-sterile mixing angle is 
small. However, the OPERA data tell us that the mixing in an overlapping L/E range cannot be small. Thus, the 
simple model we have presented does not simultaneously account for OPERA and atmospheric data. 

On the positive side, additional constraints on models of OPERA data seem to be easy to meet with our class of 
models. The coincident time of arrival of neutrinos and photons from SN 1987A is easily accommodated by making 
6*i3 sufficiently small so as to suppress v e — > v s oscillations; see Eq. (|29l) . Also, the analogue of Cherenkov radiation 
for superluminal neutrinos, pointed out in Ref. [TTj . does not apply at least to the extra-dimensional variety of our 
model since all propagation is subluminal locally, and the apparent superluminal behavior is simply a consequence 
of the bulk shortcut. (Also, with SM particles confined to the brane, there can be no Cherenkov radiation of SM 
particles from sterile neutrinos traveling in the bulk.) 

IV. CONCLUSIONS 

In light of the alleged evidence for superluminal propagation of neutrinos in OPERA data, and drawing on our 
prior work on superluminal neutrinos, we have provided herein a concrete model that accommodates superluminal 
neutrinos. We find that oscillations between active neutrinos and a sterile neutrino with a shorter geodesic path 
through extra dimensions (or equivalently, obeying a modified dispersion relation as seen from the brane) do not 
provide a simple explanation of the OPERA anomaly because of conflicts with atmospheric neutrino data. This does 
not invalidate the model in a broader context. Introducing additional sterile neutrinos and/or mixing angles into 
the framework may produce consistency with data, but we have not explored this possibility. In one variation of the 
model, the resonance may have an (Li?)-dependence Y% (instead of a simple ^-dependence) thus affording greater 
flexibility in addressing the tension with atmospheric data. It does seem, however, that an explanation of the OPERA 
anomaly using neutrino oscillations is likely to be contrived. 

Finally we remark on the generality of our results. Since the bulk-shortcut model for sterile neutrinos appears 
from the vantage point of the brane as a LIV model, our conclusions apply in generality to the larger class of models 
in which a sterile flavor state is superluminal, and transmits its greater speed to active flavor states via mixing and 
oscillations. 
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